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Abstract 

We examine in a model-independent manner the measurements that can be performed at 
B-factories with sensitivity to dark matter. If a singlet scalar, pseudo-scalar, or vector is 
present and mediates the Standard Model - dark matter interaction, it can mediate invisible 
decays of quarkonium states such as the T, J/'^ , and 77. Such scenarios have arisen in the 
context of supersymmetry, extended Higgs sectors, solutions the supersymmetric /i problem, 
and extra U{1) gauge groups from grand unified theories and string theory. Existing B- 
factories running at the T(4S') can produce lower T resonances by emitting an Initial State 
Radiation (ISR) photon. Using a combination of ISR and radiative decays, the initial 
state of an invisibly decaying quarkonium resonance can be tagged, giving sensitivity to the 
spin and CP-nature of the particle that mediates standard model-dark matter interactions. 
These measurements can discover or place strong constraints on dark matter scenarios where 
the dark matter is approximately lighter than the 6-quark. For the decay chains T{nS) — >■ 
7r'^7r"T(lS') (n=2,3) we analyze the dominant backgrounds and determine that with 400/6~^ 
collected at the T{4S), the B-factories can limit BR{T{1S) — > invisible) < 0.1%. 
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I. INTRODUCTION 



Multiple astrophysical experiments have measured the presence of Dark Matter (DM). 
The leading candidate for this DM is a particle. ^ Existing studies have concentrated on a 
single particle providing the DM density of the universe but multiple particles are allowed. 
This knowledge should lead to a systematic search for invisible decays of particles and mesons 
known to exist. However, the only particles with reported invisible branching ratios or limits 
are the vr'' and Z. 

There are at least two pieces of evidence that the DM component of the universe may be 
lighter than the Minimal Supersymmetric Standard Model (MSSM) or minimal supergrav- 
ity mediated supersymmetry breaking (mSUGRA) lightest neutralino that dominates DM 
studies. First, recent measurements of 511 keV gamma rays from the galactic center indicate 
a Gaussian profile of low-velocity positrons, Traditional MSSM or mSUGRA neutralinos 
are heavy, and their annihilation would produce too many high-energy 7-rays from neutral 
pions which decay to photons, as well as significant bremsstrahlung as their decay products 
slow until they are nearly at rest as is required explain the 511 keV line. This places such a 
scenario in strong conflict with EGRET upper limits on the higher-energy gamma flux. 
Thus, if a DM particle is responsible for the 511 keV line, it must be lighter than approxi- 
mately 100 MeV. Second, recent analysis of DM flows and caustics indicate that the CDMS 
limit and DAMA evidence for DM can be compatible due to the lower detection threshold 
of DAMA. This effect is also enhanced if there is a flow of DM through our solar system. 

Neutralinos in the general MSSM must have M^o > 6 GeV to obtain an appropriate 
relic density. 0] Constraints on even lighter DM comes from CUSB, which measured T — >■ 
7 + invisible signals, giving the best sensitivity to DM lighter than approximately 1.5 GeV, 
and losing sensitivity as the DM mass increases due to the soft spectrum of Initial State 
Radiation (ISR) photons, Modern b-factories can improve in this measurement by at least 
an order of magnitude.!^ Searches for invisibly-decaying Higgs Bosons are sensitive if the 
Higgs is heavy, has significant coupling to the Z, and the Higgs decays dominantly to dark 
matter, jl^ Finally, the LEP single-photon counting measurements limit arbitrary Standard 
Model (SM)-DM interactions; however, the Z invisible width dominates these measurements 
at LEP energies, and these experiments have no sensitivity if the SM-DM mediator does not 
couple to the electron. [lll | 

These constraints can all be avoided in many models that are not the MSSM. Two attrac- 
tive possibilities that can explain the above data are (1) a light neutralino with couplings 
to the SM mediated by a light scalar singlet. This may occur in the Next-to-Minimal 
Supersymmetric Model (NMSSM) and related models which solve the [i problem with a sin- 



glet; (2) light scalar DM coupled to the Standard Model through a new gauge boson 
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U. jisl . 14 , 0| In both cases there must be some small coupling to the Standard Model in 
order to avoid having the DM density over-close the universe, [jj It is in general possible to 
couple the DM preferentially to some quarks and/or leptons but not others. This can result 
in invisible decays of some hadrons of a given spin J and Charge x Parity (CP) eigenvalue 
but not others, and little or no signal at direct detection experiments. Building such models 
is straightforward, and several already exist in the literature. 12, 0| Our purpose in this 



letter is not to build such models, but point out several measurements that can be performed 
at colliders that are sensitive to light DM. 

At e~^e~ detectors like BaBar, Belle, and CLEO, one can use ISR to explore energy regions 
below the nominal collider energy, [ly, |l2| With one ISR photon, these experiments deliver 
the same luminosity between 9.85 GeV and 10.58 GeV as they do from interactions without 
an ISR photon at their nominal center of mass energy, 10.58 GeV. In order to identify 
the presence of a bottomonium state without observing its decay, one must require an ISR 
photon of a specific energy and/or a radiative decay. A radiative decay is any transition 
from one quarkonium state to another. We present several techniques which can be used 
to suppress backgrounds when the ISR photon is lost because it is outside the detector 
acceptance. Radiative decays from the T{4S) will have similar statistics to ISR production 
of lower T resonances, but no radiative decays of the T{AS) have yet been discovered. 

Another measurement that can be performed at B-factories is in 6 ^ s transitions such 
as + invisible and is sensitive to dark matter with masses up to 2.4 GeV, but is 

not sensifve to the J- of the mediator.Q This branehing rafo may be 50 times larger 
than is expected from the Standard Model process with neutrinos. 

The Standard Model expectation for T to invisible is r(T — > vv) = 4.14 x 10~^r(T — >■ 
e^e~) ~ 1 X 10~^ with a theoretical uncertainty of only 2 — 3%, and is sensitive to the 
bottom squark mass and R-parity violation in SUSY theories, jioj ] 

Expectations for branching ratios of hadrons into DM may be as large as a few percent. 



II. SOURCES OF DARK MATTER 



We assume that DM couples to the Standard Model through some mediating boson. On 
general model-independent grounds we expect this particle to be either a vector, scalar, or 
pseudo-scalar. If the mediator is a scalar or pseudo-scalar, only a SU (2) doublet can couple 
to hh by gauge invariance at dimension 4. This scalar doublet will generally mix with the 
Standard Model Higgs, or the CP-odd A of a Two Higgs Doublet Model. Therefore, we 
expect scalar and pseudo-scalar mediated DM to show up dominantly in interactions with 
heavy fermions such as h quarks. 

If the mediator is a vector gauge boson, giving h and h equal and opposite charge under 



3 



the new gauge group (which we assume to be a [/(I)) is sufficient to introduce the proper 



couphngs. 13|, |l5[ One need not expect this gauge boson to couple to all fermions in the 



Standard Model equally. One might expect that the first two generations are not charged 
under this new group, in order to be consistent with precise measurements of the muon and 
electron anomalous magnetic moment, as well as the lack of unexplained vector resonances 
in hadronic data. This situation would result in extremely small DM-nucleon cross sections 
for direct detection experiments. 

In order to get small DM masses in the MSSM and other models, one generally has to also 
bring down another particle mass for the purpose of getting a large enough annihilation or 
coannihilation cross section. However, in the case of DM masses less than Mx/2 = 4.73 GeV 
that we consider, one generally cannot bring down a particle charged under the SM gauge 
groups without violating existing experimental constraints.^ Having a light bottom squark 
may be one exception to this assumption; however, a recent re-analysis of available data 
indicates that this solution is now disfavored. 

The direct constraint on the annihilation mediator is the reason why existing con- 
straints 3] require MyO > 6 GeV, despite the fact that the theory is consistent with a 
massless neutralino. [22| If the mediator has some mixing with a pure singlet state, these 
constraints can be largely avoided. 

A general argument predicts that DM should be heavier than 2 GeV. j^] This is based 
on DM annihilation couplings that are proportional to Gp- However, there is nothing 
that requires DM to have something to do with weak bosons or electroweak symmetry 
breaking. The argument that DM couplings must be proportional to is based solely on 
the coincidence that the DM annihilation cross section (c.f. Eq. |3J is similar in size to weak 
cross sections. This may simply be a numeric coincidence and annihilation cross sections 
need not be proportional to Gp- If we simply assume that DM exists and < 6 GeV is 
an allowed region, we are forced to recognize that it must have picobarn cross sections with 
some Standard Model particle. These expected cross sections are explored in the following 
section. 



III. DARK MATTER COUPLING EXPECTATIONS 



The required rate of DM annihilation can be naively estimated. We will not accurately 
compute the relic density since we are not proposing a specific DM model, but one can get 



The usual particles that are made light in supersymmetric models are the stau or a Higgs. The stau is 
often the next-to-lightest supersymmetric particle and undergoes t-channel coannihilation with the LSP. 
The Higgs mediates s-channel annihilation when there is significant higgsino or wino fraction in the LSP. 
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an order of magnitude estimate for s-wave annihilation using 



2 O.lpb-c 

Qxh - — — — . (1) 

{av) 

Where Qx = Px/Pc is the rehc density for species X relative to the critical density pc, h is 
the Hubble constant, and {av) is the thermally averaged annihilation cross section of the 
DM into Standard Model particles. Using the central value of the WMAP result for 
Qxh"^ = 0.113, we can invert this equation and solve for the required annihilation cross 
section for light relics 

(av) = 0.88pb. (2) 

The velocity v appearing here is the M0ller velocity, the relative velocity of annihilating 
particles at the temperature they froze-out. The approximate temperature at freeze-out is 
T = m^/xpo where is the mass of the DM and xpo is an expansion parameter evaluated 
at the freeze-out temperature that is Xi^o ~ 20 — 25 depending on the model. Thus the 
average momentum for a fermion is /c^T and therefore the average relative velocity is roughly 
l/xpo- Foi' a; = 20 at freeze-out we have: 

a{xx ^ SM) ^ 18pb. (3) 

The invisible branching ratio of a hadron can then be estimated by assuming that the 
time-reversed reaction is the same, a{ff — ^ xx) — ^ixx ~^ ff)- This assumption holds if 
c^rrtf and Mx — 4m^+6m^Tpo- We assume that the DM mediator is not flavor changing 
and that annihilation occurs in the s channel. ^ Therefore, the best-motivated hadrons to 
have an invisible width are same-flavor quark-antiquark bound states (quarkonia). The 
CERN Yellow Report provides a thorough review of quarkonium physics, 

The invisible width of a hadron composed dominantly of qq is given approximately by: 

V{H ^ XX) = flMHa{qq ^ xx) (4) 

where fn is the hadronic form factor for the state H, and Mh is the hadron's mass. 

We can predict an approximate expectation for the branching ratios for narrow states. 
Some of the most promising are: 

BR{T{IS) xx) ^ 0.41% BR{J/'^ xx) ^ 0.023% BR{ri xx) ^ 0.033% 

(5) 

Branching ratios for scalars and pseudo-scalars tend to be smaller since those states are 
wider. This estimate does not take into account kinematic factors arising from the mediator 
mass and DM mass. These factors can both enhance or suppress these branching ratios. 



^ A i-channcl mediator is possible, but this requires that the mediator carry color and electromagnetic 
charge, and therefore is unlikely if we consider < 5 GeV. 
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If a particular hadron H decays invisibly, then that hadron must mix into the mediator 
M before decaying into Dark Matter. If M does not violate the discrete symmetries C 
and P, H and M mix only if they share the same spin, C, and P eigenvalues. Therefore, 
the observation of an invisible decay not only constrains the mass of the dark matter and 
mediator, but may also uniquely identify the spin, C, and P of the mediator. An invisible 
decay does not have sensitivity to the spin of the Dark Matter itself. 

With running B-factories BaBar and Belle having roughly 400 fb~^ recorded, these ex- 
periments may already have tens of thousands of DM production events, if the DM is 
kinematically accessible. 

IV. BOTTOMONIUM PRODUCTION VIA ISR 

Bottomonia can be identified by observing the particles emitted when it makes radiative 
transitions to lighter bottomonia. These transitions are show in Fig. ^ Since the states are 
fairly narrow, the energy of the photon radiated in a transition (or kinematics of a particle 
pair) gives a clean way to select specific quarkonia transitions. The CLEO experiment has 
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FIG. 1: spectra for bottomonium, with the spin state on the horizontal axis 
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provided measurements of most of the quarkonia transitions. j25| The number of ISR T 
production events collected by the BaBar and Belle experiments is now competitive with 
that collected by the traditional method of scanning the resonance. Furthermore, off-peak 
data can be used for the measurements we propose. 

The ISR cross section for a particular final state /, with e^e~ cross section is to 

first order 

da{s, x) 



dx 



W{s,x)-af{s{l-x)) (6) 



where x = Ej is the energy of the ISR photon in the CM frame, and ^/s is the CM 
energy. The function 



W{s,x) = (3 



(7) 



describes the energy spectrum of the ISR photons, where (3 = ^(2 In ^ — 1) and S take into 
account vertex and self-energy corrections; a is the electromagnetic coupling constant and 
rrie is the mass of the electron. At the T(4S') energy, f] = 0.088 and 6 = 0.067. By tagging 
the ISR photon, B-factories can explore all of the vector bb bound states. 

This function (Eq. [7j) is strongly peaked in the forward and backward directions, so ISR 
photons will be close to the beamline. Assuming a detector acceptance of —0.9 < cos(^) < 
0.9, the fraction of the nominal luminosity at the T(45') resonance delivered to the T(IS'), 
T(2^) and T(3^) resonances is 1.9 x 10"^ 3.2 x 10^^ and 5.0 x 10^^ respectively. This 
results in hundreds of thousands of events per resonance with current recorded luminosities. 

If one does not require that the ISR photon is identified because it is not in the detector 
acceptance, the production of lower resonances is larger by a factor 5 to 7. Q| The fraction of 
the nominal luminosity at the T(4S') delivered to the T(IS'), T(25'), and T(3S') resonances 
is 8.5 X 10-^ 1.5 X 10-^ and 2.3 x 10"^ 

One can then tag quarkonium states by looking for a particular radiative transition such 
as T(25') Xbo{^P)l or T(35') T(15')7r+7r^. One can examine the decay modes of the 
tagged state in a decay-mode independent manner. 



V. THE KINEMATICS OF ISR PRODUCTION WITH RADIATIVE DECAYS 

The irreducible physics background for these invisible decays coming from a pair of neu- 
trinos and pions is extremely small due to the fact that weak cross sections are suppressed 
by (Mx/Miy)^ ^ 0.01, the final state has high multiplicity, and our signal has resonant 
enhancement. For example a{e^e^ ir'^n^uu) ^ 10^^ pb before applying any cuts. There- 
fore, the dominant backgrounds will come from unrelated processes that do not actually 
have a neutrino or DM, or neutrinos from r decays. 
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The knowledge that resonances are formed in our signal gives us the kinematic constraint 
that the square of the four-momenta forming a resonance must be equal to the resonance 
mass-squared. For the production of a single T resonance via ISR and not observing the 
ISR photon, there are 2 kinematic variables that are undetermined. These variables are 
associated with the four- vectors of the ISR photon and T which we presume decays invisibly. 
With n intermediate resonances decaying radiatively to each other, in the event that the 
ISR photon is unobserved, we can predict all but 2 — noi these undetermined variables. The 
most important radiative decays of bottomonium are listed in Appendix sorted by cross 
section. Up to two intermediate resonances can be created with sizeable event rates. 



A. Resonance Constraints 



We can either observe the ISR photon, or require that its angle with respect to the 
beamline is consistent with it being outside the detector acceptance using the kinematic 
constraints. 

With one intermediate resonance and one radiative decay (e.g. T(2S) T(15')7r+7r~), 
we define Mi (e.g Mx{25)) to be the mass of the intermediate resonance and M2 (e.g. Mx(i5)) 
to be the mass of the final invisibly-decaying state. We can predict all but one kinematic 
variable using the measurement of the particles emitted in the radiative transition and the 
beam constraint. In the center-of-mass frame: 

where = {Er]Pr) is the sum of the four momenta of all the particles emitted in the 
radiative transition, = Pr^iPr, Pr = \Pr\, and s = M^f^^g--^ is the center-of-mass energy. 
Our predicted angle 6 is the angle between the ISR photon and pr- cos 6* can be related to 
AM^ = Ml — M|. However, when expressed as an angle it is clear that it can still be used 
when the ISR photon is unobserved. 

We can also invert Eq. |S1 to come up with a cut on the energy of the radiated system 

which is useful for single photon transitions. 

Backgrounds without resonances can be distributed outside the physical region — 1.0 < 
cos 6* < 1.0 since for background without resonances, Eq. |Hldoes not describe any physical 
angle at all. Therefore, a cut requiring —1.0 < cosO < 1.0 will suppress most backgrounds 
by a factor 10^ — 10^ (c.f. Table. H}. If the ISR photon is unobserved we also know that it 
must be outside the detector's acceptance. Assuming the EM calorimeter extends between 
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^min and Omax as 86611 from th6 c6nt6r-of-mass-fram6, th6 angl6 b6tw66n th6 ISR photon and 
th6 b6amlin6 (6'isr) must satisfy ^isr < 9min or 9isr > Omax- Since the angle with r6sp6ct to 
th6 b6amhn6 6r of Pr is m6asur6d, this amounts to th6 r6striction: 

\9 — 9r\ < 9min OI \9 + 9r\ > 9„iax (10) 

Furthermor6, th6 signal p6aks in both th6S6 r6gions as {9 ± 9^)^ . This corr6sponds to th6 
ISR photon b6ing n6arly parall6l to on6 of th6 b6amhnes. 

If th6 ISR photon is observed both -Eisr and cos^ ar6 measured and can be directly 
compared to (jH)). The final kinematic variable can be taken to be the angle between the 
plane defined by the beamline and pr, and the plane defined by beamline and pisR- Only 
if the ISR photon is observed can this angle be determined. This angle will only provide 
power in suppressing background if the background happens to peak in this variable. If the 
ISR photon is unobserved, this angle is not knowable in principle. 

When two intermediate resonances are formed, and two radiative decays occur, all kine- 
matic variables can be determined by measuring the energy and momentum of the particles 
in the radiative decay. The constraints above still apply for the first transition. The new 
constraint available allows us to predict the angle between the ISR photon and the second 
radiative decay: 

v/i M|-M| + r,^ + 2ri-r2 s + 

COSt^ = -— TT ^ 77^ 777 (H) 

We can further apply the same trick as with cos 9 to require the ISR photon to be in the 
beam line, replacing 9 with 9' in Eq. ^| 

It should be noted that it is possible to emit two or more ISR photons. In this case, 
Eqs. IHlim are accurate in the limit that the ISR photons are coUinear. 



B. Missing Momentum and Standard Model Decay Backgrounds 

Considering that there are invisible particles in the final states we are interested in, there 
is an irreducible background when the final state decays to visible Standard Model particles, 
but those particles lie outside the detector acceptance. A missing momentum cut can force 
the transverse momentum of the final state to be large enough that one can ensure that 
its decay products would be in the detector volume. However, this is almost completely 
useless for machines running at the T{4S) for the following reason: assume the final state 
undergoes a 2-body decay, those decay products lie exactly on the edge of a symmetric 
CLEO-like detector at 6'min and 6'niax = tt — ^min, and the final state is at rest at the center. 
The cut required on the transverse component of the sum of radiative transition particles is: 

{Pr)T > \{EcM - EiSR - Er) siu 9^in ^ 2GeV, (12) 
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which is larger than the visible energy {Ejsr + Er) in any radiative transition. Here Eqm 
is the center of mass energy, -Eisr is the energy of the ISR photon, and E^. is the sum of 
energies of particles emitted in the transition(s). This background is further discussed in 
Sec. Ivm 

For a collider running at O{30GeV), requiring the ISR photon to be visible provides 
enough transverse momentum that the decay products of the final state must lie in the 
detector acceptance. However, the ISR production cross section of T's is reduced by a 
factor ~ 100. This requirement also eliminates most of the two-photon background, as 
discussed in Sec. IVI Al 



VI. INVISIBLE UPSILON DECAYS 



To demonstrate that invisible widths can be measured using ISR and radiative decays, we 
concentrate on the modes T{nS) T(lS')7r+7r~ (n=2,3) for colliders running at the T{4S) 
since these modes have the largest cross section (2.91pb, 0.784 pb). Many decays^ have the 
visible topology vr+vr" + wj, n > 0, which might be useful for triggering. These modes seem 
the most promising since the running B-factories BaBar and Belle do not have the sensitiv- 
ity in their calorimeter that CLEO had, and they have excellent c harg ed particle tracking 
resolution by design. These transitions were measured by CLE0.[l6| A table of possible 
transitions and their tagging signatures is presented in Appendix ^ The T{3S) mode has 
harder pions, and therefore may have a higher reconstruction or triggering efficiency than 
the T(2^) mode. 

For background simulations we have used PYTHIA j3l and CompHEP with r decays 
simulated with TAUOLA 28[. We have smeared charged tracks according to the BaBar 
detector resolution 

^ = 0.21% © lA%pt, (13) 
Pt 

charged tracks must have pt > 100 MeV, and all objects must lie within the detector 
—0.87 < cos^ < 0.96 from the center-of-mass frame (i.e. BaBar geometry). For photons we 
require E > 20 MeV and smear their energy according to 

(Te 1.0% 3mr 

^ = 1.2%©— -jj, ag = a^ = 2mr© , 14 

E {E/GeVf^ ^E/GeV 

Finally, for charged tracks we do not differentiate 7r+, e~^, /i"*" or and assign each charged 
track to have a mass after smearing its momentum, since tracking information is reliable 



3 For instance T{3S) ^ T(2S')77 ^ T(lS')777r+7r-, T(3S') ^ T(2S') ^ Tr°n° i:{lS)n°TT°TT+n^, 
T{SS) ^ /ib(lP)^+7r- ^ 77fc(15)^+^-7, T(3^) ^ X6o(2P)7 ^ Vb{lShr^, T{3S) ^ T(2^)^+^- ^ 
X&o(l-P)7r+7r-7. 
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but particle ID is not. These tracks generally are soft enough that they do not enter the 
calorimeter, or enter at a grazing angle. 

The cuts proposed in Sec. IV Al are drastic and in general can result in a background 
suppression of 10^ or more, including effects of detector resolution. Due to the large size of 
these backgrounds, detailed detector resolution effects and multiple scattering will be very 
important. Therefore, each background should be measured directly in order to estimate 
the signal contamination. Therefore, we present estimates of each background and the 
level to which they can be suppressed including smearing. The exact numbers may change 
significantly when detector effects are taken into account. 

Aside from the purely kinematic constraints presented in Sec. IV A| there is no further 
angular information from the matrix element that is useful to identify the signal. In the case 
of a final state T from a di-vr transition, there is no spin correlation between the outgoing 
pions and the final T. The operator involving a coupling between the polarization of the 
T and the momenta of the pions is D-wave suppressed and measured to be very small, [soj ] 
Furthermore, the vrvr invariant mass spectrum in the T(3S') T{lS)7nr transition is not 
well understood theoretically and may involve another intermediate state. ^ 

Including all backgrounds and a realistic smearing of energy and momenta, the measure- 
ments proposed can limit at 2a sensitivity 

5i?(T(15) ^ invisible) < 0.113% (15) 

using the T{2S) mode and 

BR{T{1S) invisible) < 0.335% (16) 

using the T{3S) mode. The combined 2a sensitivity is then 

BR{T{1S) ^ invisible) < 0.107%. (17) 

The dominant backgrounds are discussed in the following subsections. Their numeric 
importance and cuts needed to suppress them is summarized in Table H] 



A. Photon Fusion Background 



The two photon fusion process occurs when both incoming beams emit a photon and 
those photons annihilate into electrons, muons, taus, or hadrons. This cross section is very 
large, in the hundreds of nanobarns. Furthermore, our signal spans the region < < 1 
GeV^ in which non-perturbative QCD effects dominate hadron production. Due to this, a 
reliable simulation of hadron production is not possible and in any case should not be relied 
upon due to non-perturbative effects. This background must be measured directly. 



11 



Backgrounds to T{3S) T(lS')7r"'"7r 
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TABLE I: Dominant backgrounds to the processes e"'"e~ T{nS)jjsR T{IS)'it~^tt~jjsr- Cuts 
are shown in the left column, followed by the effective remaining cross section in each background 
channel after the cut. Each line includes all the cuts above it. The photon is assumed visible, 
and any charged track is assumed to be a pion. For the T(25) mode, the 2a sensitivity on 
BR(T(1S') ^ invisible) ~ 0.113%; for the T(35) mode, the 2a sensitivity is ~ 0.335%. 

To demonstrate that this background can be overcome, we simulate 10 fb~^ of the lepton^ 
production processes e~^e~ — > e^e~l~^l~ , e^e~ — e^e~ + hadrons and e~^e~ — e^e~l~^l~^. 
We simulate the first two backgrounds using PYTHIA, and the second using CompHEP ^ in 
the equivalent photon approximation, [s^l At these low energies, ir/fi separation is generally 
unreliable since the muon is not energetic enough to reach the outer muon detector. The 
l^l^ cross section is also about an order of magnitude larger than the vr+vr" cross section, 
making the the most important background in any case. The fi'^fi" cross section is 38.5 
nb for the BaBar detector geometry, assuming both muons are visible in the detector volume. 
This is sufficiently large that it overwhelms the physics signal that is normally triggered on 
at the B-factories (about 1 nb). Therefore, the rate must be reduced at the trigger level. 
Requiring an extra visible photon reduces this cross section to a triggerable level. In the 
case of the e+e~ — > e^e~l^l~ background, an extra photon comes from Initial/Final State 
Radiation. 

Di-lepton events produced in photon fusion have the characteristic that the leptons are 
back-to-back in the plane perpendicular to the beamline. By contrast, our signal is a 3- 
body decay, so only a small fraction are back-to-back. Therefore, = and cos 6' (c.f. 
Eq. IHI) will be very large. As can be seen in Table IJ a cut on cos 9 alone can remove this 
background. If detector resolution effects cause this background to bleed into the signal 
region, a cut < vr, where is the angle between the leptons in the plane perpendicular 

^ Here lepton refers to an e or /i. 

^ The 77 process is not included in PYTHIA. 
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to the beamline can also remove this background. Therefore, this background should be 
carefully studied. This can be done by identifying singly-tagged photon fusion events, where 
one of the initial electrons is deflected into the detector. Experiments such as DELPHI have 
also employed far-forward particle detectors to identify the photon fusion signal when the 
electrons are deflected by a small angle (known as the Small Angle Tagger and Very Small 
Angle Tagger). 

B. di-r Background 

The only irreducible physics background to these processes that have true neutrinos comes 
from r decays. The dominant source is e^e~ t~^t~ via a virtual photon, where both r's 
decay to pions. The total t~^t~ cross section is 993 pb at tree level, r's can also be produced 
in photon fusion, with a cross-section of approximately 22 pb. 

For example, we consider the di-r background to the process T(3S') T(15')7r"''7r~. We 
require exactly two charged tracks and one photon visible in the detector. Generally the 
photon comes from a ir^ decay in which the other photon is outside the detector's acceptance, 
or final/initial state radiation. 

The effect of the kinematic cuts proposed in Sec. IV Al are shown in Table |lj The di-r 
background generally has very different kinematics than our signal, as well as extra vr'^'s. It 
can be reduced below 0.1 fb with these cuts. 

C. Two-Body Decay Background 

A background to all processes is true resonance production where the final state resonance 
decays via any 2-body decay and its decay products lie outside the detector acceptance. This 
background is irreducible, but is accurately measured in events with both the radiative decay 
of interest and the final hadron decaying to a 2-body state. This gives roughly 10 times the 
statistics on measuring this background, so it can be subtracted. 

For 2-body decays this amounts to an irreducible background that has a branching ratio 
f2^l. Here /2 is the fraction of 2-body final-state bottomonium decays plus multi-body 
decays which are arranged such that all decay products are outside the detector acceptance. 
Q is the fraction of the solid angle covered by the detector. This background can be directly 
measured by relying on the sample of non-invisibly decaying final state particles provided 
by the ISR + radiative transition technique. We take /2 = 5% and Q = 91.5% for the 
BaBar detector geometry. The final decay is uncorrelated to its production mechanism 
and the radiative transition, and therefore is isotropic in detector angle 6 and 0. Events 
with a radiative transition and visible bottomonium decay give a measurement of all decay 
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channels of the final state particle, including effects of detector resolution for the radiative 
transition, ISR smear, and multi-body decays. It should be noted that one cannot simply 
take this background sample and pretend the beamline bisects the detector in a different 
direction. The asymmetric boost of modern B-factories changes the size and area of the 
would-be beamline and this must be taken into account. 

D. Drell-Yan 

Direct production of vv is small since the neutrinos must come from a Z or ^ which 
are heavy. For instance, BR(T(1S') vv) ~ 1 x 10~^. 0] Only the vector resonances have 
a sizeable branching fraction to neutrinos, since they can mix directly with the Z . Scalars 
and pseudo-scalars can only emit neutrinos in loop suppressed processes. 

Modern B-factories do not have the sensitivity to test this branching ratio, and therefore 
it is not a background. 

VII. CONCLUSIONS 

Measurements of invisible branching ratios of mesons are extremely important, given the 
established evidence for Dark Matter (DM) and the knowledge that most DM scenarios 
require some Standard Model-DM interaction. Given tight constraints on flavor changing 
neutral currents, the most important mesons to examine are flavor neutral bound states 
of quarks. Running high-luminosity B-factories motivate looking for invisible decays of 
bottomonium flrst. 

ISR and radiative decays provide a powerful method to measure the invisible branching 
fractions of the bottomonium resonances. If the DM is lighter than Mx/2, annihilation 
of DM into standard model particles is expected to have a picobarn-scale cross section. 
While the sensitivity achievable is not capable of measuring the Standard Model T — > vv^ 
decays to dark matter should be signiflcantly stronger than decays to neutrinos, due to 
the {My/MzY suppression of the Standard Model process. These techniques can limit 
BR{T{1S) — > invisible) < 0.1%, which is sensitive enough to discover dark matter if it 
couples in this manner. 

DM with a mass < 5 GeV is generally allowed in models. Direct detection experi- 
ments are very insensitive in this mass region, and would also be insensitive if dark matter 
preferentially couples to heavy quarks. Therefore, alternative methods to discover DM are 
required if the DM is this light. We strongly encourage experimental teams at BaBar, Belle, 
and CLEO to pursue these techniques. 
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APPENDIX A: BOTTOMONIUM EVENT RATES 



In the following we give the expected cross section for bottomonium production assuming 
Eqm = ^T(45) = 10.58 GeV. It should be noted that both on-peak and off-peak data can 
be used for this analysis. 

For each final state, a "tagging topology" is given, which is the set of particles visible 
in the detector's acceptance. In each case the particles in the tagging topology have a 
well-defined kinematics as outlined in Sec. IV Al In cases where 7isr is not listed in the 
tagging topology, the cross section includes both when the ISR photon is visible, and when 
it lies outside the detector acceptance. When 7isr is listed, the cross section corresponds to 
requiring the ISR photon to be visible for the BaBar detector geometry, —0.87 < cos 6 < 0.96 
in the center-of-mass frame. 

We catalog only existing, measured resonances and transitions of bb quarkonia here (with 
the exception of the undiscovered r]b). Other decay chains will be possible involving radiative 
decays of the T{4S) when those decays are discovered. This clean method of tagging the 
initial state will allow the discovery and cataloging of many more radiative decays, improving 
statistics from what is listed below. 
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A-1. Vector Mediated Dark Matter 



Final state 



Decay chain 



Tagging topology 



T{1S) 



T{2S) T{1S) 
T{2S) T{1S) 
T{3S) T{1S) 
T(25)^Xw(lP)^T(lS) 

T(35) ^ T{1S) 
T(25)^X62(1P)^T(15) 
T{3S)^Xbi{2P) 
T{3S) T{2S) 
T(35) ^ Xb2{2P) 
T{3S) T{2S) 
T(35) ^ T(25) 
T(35) ^ T(25) 
T(35) ^ T(25) ^ T(15) 
T(3^)^Xw(2P)^T(15) 
T(3^) T{2S) 
T{3S) ^ X62(2P) 
T(35) xw(2P) 



^ T(15) 
T(15) 
^ T(15) 
T(15) 
T(15) 
T(15) 



T{1S) 

> T(15) 

> T(15) 



7lSR 

7 

7 
7 

77 
7 

77 



7 



7 

7 
7 

7 



TT ' TT 



7 
7 



7 



Total 



T(2S) 



T(35) ^ T(25) 
T(3^) ^ T(2^) 
T(35)^Xw(2P)^T(25) 

T(35) ^ T(25) 
T(35)^Xfe2(2P)-^T(2^) 
T(35)^X6o(2P)^T(25) 



7lSR 

77 

7 

7 
7 



7 

7 
7 



Total 



A- 2. Pseudoscalar Mediated Dark Matter 



Final state 



Decay chain 



Tagging topology 



^6(1^) 



Total 



T(3^) ^ h{lP) ^ VbilS) 
T{3S)^hilP)^VbilS) 
T(35) ^ X6o(2P) ^ r/fe(15) 





7 




7 


7 
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A-3. Scalar Mediated Dark Matter 



Final state 


Decay chain 


Tagging topology 


o-T(pb) 


Xw{lP) 


T(25) ^ XboilP) 
T{3S) ^ T{2S) ^ xmilP) 
T{3S) ^ T{2S) ^ XboilP) 
T{3S) ^ T{2S) ^ XboilP) 


7iSR 7 

77 7 

Tr^vr*^ 7 


0.0937 
0.0333 
0.0186 
0.0133 


Total 




0.159 


Xw(2P) 


T(35) ^ X6o(2P) 


7iSR 7 


0.188 


Total 




0.188 
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